Prediction and Measurement of Transient Responses of First Difference Based Chaos 
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Chaotic behavior can be produced from difference equations with unstable fixed points. Differ- 
ence equations can be used for algorithms to control the chaotic behavior by perturbing a system 
parameter using feedback based on the first difference of the system value. This results in a system 
of nonlinear first order difference equations whose stable fixed point is the controlled chaotic be- 
havior. Basing the feedback on the first difference produces distinctly different transient responses 
than when basing feedback on the error from the fixed point. Analog electronic circuits provide the 
experimental system for testing the chaos control algorithm. The circuits are low-cost, relatively 
easy to construct, and therefore provide a useful transition towards more specialized real- world ap- 
plications. Here we present predictions and experimental results for the transient responses of a first 
difference based feedback control method applied to a chaotic finite difference 1-dimensional map. 
The experimental results are in good agreement with predictions, showing a variety of behaviors for 
the transient response, including erratic appearing non-steady convergence. 
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I. INTRODUCTION 

May showed that the logistic equation, a recursion re- 
lation based on a quadratic return map, can generate 
chaotic behavior. 1 Proportional feedback methods can 
control chaos in these systems^ - — Typically in these con- 
trol methods a system parameter is perturbed by an 
amount proportional to the "error", the difference be- 
tween the current system value and the unstable fixed 
point. This control method is simple proportional feed- 
back (SPF). Pyragas introduced an alternative in which 
the perturbation is proportional to the first difference of 
the system value (difference between current and previ- 
ous values), and therefore does not depend on the un- 
stable fixed points This delayed feedback control (DFC) 
method successfully controls chaotic behavior in a variety 
of cases.— Here we apply DFC to a chaotic finite differ- 
ence return map and derive the mathematical form of the 
transient response and verify it experimentally with an 
analog electronic circuit. 

The system whose chaotic behavior is to be stabilized 
is a 1-dimensional finite difference map for system value 
x with system parameter Oq 



t-n+l 



(1) 



Stabilization is accomplished using DFC by perturbing 
parameter ao by an amount proportional to the first dif- 
ference, d(x n ) = (x„ — x n _i). Thus when stabilization 
is turned on we have a system of two first order nonlinear 
finite difference equations 



/ (^*roj a n ) 

a n +\ = a + K (/ (x n ,a n ) 



(2a) 
(2b) 



K is the gain for the feedback that perturbs system pa- 
rameter a n from ao- The first difference can be written 
in terms of the error term 



d(x T , 



_A X vy —A X f. 



(3) 



where the error is Ax = x — x* and the fixed point is x* . 
Thus DFC has the desired property of the perturbation 
vanishing when the system value attains the fixed value. 
Stabilization is attempted only when the magnitude of 
d{x n ) is within a specified window for control. 

We perform a linear stability analysis to determine the 
conditions for stabilization of the chaotic behavior and to 
calculate the form of the transient response. The control 
algorithm is implemented in an analog electronic circuit. 
Measurements from the circuit compare well with pre- 
dictions. A variety of transient responses are found de- 
pending on the gain of the feedback. Some show a steady 
convergence of the system value to the fixed value, while 
other cases appear to converge somewhat erratically. The 
behaviors are distinctly different from the transient re- 
sponses obtained when using SPF control based on the 
error Ax. The electronic circuit provides a real- world 
system for testing implementation of the control algo- 
rithm. 



II. PREDICTIONS 

Here we do linear stability analysis in order to deter- 
mine the conditions under which chaotic behavior can be 
stabilized and we find the form of the transient response 
when stabilization is turned on. We use the first differ- 
ence d(x n ) as our variable of interest since the perturba- 
tion of system parameter a n is Aa„ = Kd(x n ). Using 
d(x n ) can be useful in real- world applications since the 
system values are easily obtainable, whereas the fixed 
value may be unknown or may drift. 

The fixed point of Eq. @ satisfies x* = f(x*,a ) and 
a* = ao. The Jacobian for linearization of Eq. ([2]) about 
the fixed point is 

j =U(Z"-d4) (4 » 

where f x and f a are the partial derivatives evaluated at 
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the fixed point. We point out that knowledge of the 
fixed point is required to evaluate the partial derivatives. 
However during the attempt to control chaotic behavior 
it is not necessary to monitor and update a potentially 
drifting fixed point, as long as the partial derivatives do 
not change too much. In contrast, a feedback method 
based on the difference from the system value to the fixed 
point requires monitoring and updating if the system's 
fixed point can drift. 

The trace and determinant of the Jacobian are (f x + 
Kf a ) and Kf a , respectively. These give the condition for 
a stable fixed point, 



where 



\fx + Kf a \ <l + Kf a <2. 



(5) 



The right-side inequality requires Kf a < 1. We are in- 
terested in situations when the fixed point is unstable, so 
f x < —1. Therefore Kf a > 0, and the left-side condition 
requires {—f x — Kf a ) < (1 + Kf a ). Thus the condition 
on K f a for a stable fixed point is 



7x-i 



< Kf a < 1. 



(6) 



This gives the result that DFC can not control chaos if 
the slope of the return map at the fixed point is steeper 
than —3. 

Now we find the form of the transient response of the 
first difference. It follows from Eq. ^ that Eq. Q is the 
Jacobian for d(x n ). Let y n = d(x n ) and look for solutions 
y n = A™. Characteristic multipliers are 



, fx + Kfa 

A± = ± 



'fx+Kfa 



giving solution 



Vn = ciA; + C 2 A" 



- Kfa (7) 



(8) 



The coefficients c\ and C2 are determined by two consec- 
utive observed values yo = c\ +C2 and y% = c\\ + + C2A_. 
Solving for the coefficients gives 



Cl 



'■'2 = 



-X-yp + yi 
A+ - A_ 
A+j/o - y\ 



A_ 



(9a) 
(9b) 



Setting the discriminant in Eq. to zero gives the 
value of K fa for transition from geometric to sinusoidal 
behavior, 



Kf a =2- f x - 2y/l - f x 



(10) 



For values of Kf a larger than in Eq. (|10[) and less than 
one the characteristic multipliers are complex conjugates 
and the convergence is sinusoidal. Using the Euler rep- 
resentation, A± = re ±l< ^, and taking the form of Eq. (j8]) 
the solution is 



y„ 



C [1 



(re 



-i<t>\ 



(11) 



yo , . ( yoRe(x+) - yi 



+ i 



2 V 2iro(A+) 
Thus the solution for the first difference is 



c e 



w 



(12) 



Vn = c e lS (re^)" + c e- ie (re"^)" = 2c r™ cos [n<f> + 6) 

(13) 

where 



VWa, 



= arctan 



( 4Kf a - (fx + Kfa) 



and 



9 = arctan 



fx + Kfa 



y Re(X + ) - yi 
y Im(X+) 



(14) 



(15) 



(16) 



The cosine term can cause erratic appearing convergence 
distinctly different from the steady geometric conver- 
gence for the real multipliers in Eq. 

The result is that we predict convergence (control) 
for feedback gain K between (1 + f x ) / (—2f a ) and l// a , 
with Eq. (1101) giving the transition from geometric con- 
vergence [Eq. to sinusoidal convergerence [Eq. (IT3|) ]. 

As an example we consider the 1-dimensional Henon 
map: 



with fixed point 



/ (x, a) = 1 — ax 2 



-1 + Vl + 4a 
2a 



and partial derivatives evaluated at (x*,ao) 
fx = -2a x* = 1 - \f\ + 4a 



fa 



( x *y 



(17) 



(18) 



(19) 



(20) 



For values of ao between 1.4 and 2 Eq. <JTTJ) displays a 
variety of unstable behavior including high period oscil- 
lations and chaos J£ For ao = 1.9 we find f x = —1.93 and 
f a = —0.259, so convergence of d(x n ) is predicted for val- 
ues of K between -1.8 and -3.87 with the transition from 
geometric to sinusoidal convergence at -1.96. 



III. CIRCUIT AND MEASUREMENTS 

Figure [T] shows the circuitry used to apply the control 
algorithm to a function block circuit / (x, a) that per- 
forms analog computation of a chaotic return map. The 
voltage V corresponds to the system value x, where the 
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FIG. 1: The circuit for controlling chaotic behavior of the return map x„+i = f(x n ,a n ). Op amps are LF412. 



scaling factor of the AD633 multiplier integrated circuit 
must be used so that x n — V n /(IQ volts). Here we use the 
function block circuit shown in Fig. [5] that calculates the 
1-dimcnsional Henon map Eq. (fT7)) ^ We have also used 
function block circuits that produce the Logistic map and 
the tent map. 

At the upper left in Fig. Q] the unperturbed param- 
eter value ao is added to the perturbation Aa„ to cre- 
ate system parameter a n . This is input, along with sys- 




FIG. 2: Henon function block circuit for f (x,a). Relation 
between voltage and system value is x — V/10. The xlO 
noninverting amplifier at the output accounts for both a and 
the +1 not being multiplied by 10, the scaling factor of the 
AD633 multiplier. 



tern value voltage V n , to the / (x, a) circuit block which 
produces the next system value voltage V n +\. The sub- 
traction op amp creates the first difference AV^+i = 
V n jf-i — V n that is used to create the perturbation for 
the next iteration, Ao n +i- A^ n+ i is passed to an abso- 
lute value/comparator stage and to a gain stage which 
produces Aa n+ i. The output of the comparison stage 
(LM339) controls the gate of the FET in the gain stage 
so that if | Ay | is larger than the control window then the 
gate goes to —5 volts turning off the FET and thereby 
setting feedback gain K = 0. A nonzero value for K 
is determined by the inverting op amp adjacent to the 
FET. For the values shown, 47fcfi and 13M3, K = -3.6. 
The sign of K is easily switched by changing the order of 
inputs V n and V n +i to the subtraction amplifier. Prior 
to the FET |AV| is divided by 10, the scaling factor of 
the AD633 multiplier used in the / (x, a) circuit block, 
to convert from voltage |AV| to |Ax|. The sample/holds 
(LF398) act as shift-registers for the iteration of n un- 
der the control of the 555 timer circuit. With the 68fcil 
and 0.001/i/ shown in the schematic the period is about 
100/is. 

Data were collected with a Tektronix TDS 3000 oscillo- 
scope. The control circuit was periodically gated on and 
off (circuitry not shown) so that it was possible to trig- 
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ger from the gating signal in order to capture the entire 
control of chaotic behavior. Figures [3] and [U show data 
and the effect of the gating. Also apparent is the control 
window. Control was gated on at t = in both cases, 
but in Fig. [3] |AV| was not within the control window 
until about t — 4.5 ms. 




--8.005 0.000 0.005 0.010 0.015 0.020 

Time (sec) 

FIG. 3: Data for control of chaotic Henon system with feed- 
back gain K = —1.96, a value predicted to result in steady 
convergence of the system value. Control was gated on at 
t = and off at t = 0.012. Also shown is the paramater value 
a — do + Aa with ao = 1.9. 



IV. RESULTS AND DISCUSSION 

Figures |3] and @] show the effectiveness of the control 
circuit in Fig. [T] when applied to the Henon circuit with 
parameter ao = 1.9, a value that gives chaotic behavior. 
The figures show the measured voltages for the system 
values V n = 10x n and parameter values a n . Figure [3] uses 
feedback gain K = —1.96 corresponding to the transition 
between geometric and sinusoidal convergence [Eq. (I10[) ]. 
and Fig. 2] uses K — —3.53 corresponding to non-steady 
sinusoidal convergence. 

The DFC method uses feedback proportional to the 
first difference of system values, d(x n ) = x n — x n —\- 
When the system is successfully controlled d(x n ) — s- 0. 
In Section II we showed that the nature of the conver- 
gence of d(x n ) depends on the feedback gain K. The ap- 
proach to zero for d(x n ) may be steady [Eq. p])] or may 
appear somewhat erratic [Eq. dTBl ], Figures [SJ [BJ and [7] 
show data and prediction for y n (actually AV n — 10y n ) 
for three values of K showing the variety of convergence. 
Predictions were made by using two successive measured 
system value differences for yo and y\ in Eqs. © and 
to determine the coefficients for Eqs. © or (|13l) . 
The convergence in Fig. [5] is steady, in Fig. [6] it appears 
somewhat erratic, and in Fig. [7] a pattern is apparent al- 
though the convergence is not steady. In all cases there 
is good agreement between the prediction and measure- 
ment. 
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FIG. 4: Data for control of chaotic Henon system with feed- 
back gain K — —3.53, a value predicted to result in non- 
steady convergence of the system value. Control was gated 
on at t = and off at t — 0.012. Also shown is the paramater 
value a — ao + Aa with ao = 1.9. 
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FIG. 5: Data (dots) and prediction (open circles) showing the 
first difference of system values d(xi) = (xi —Xi—i) for steady 
geometric convergence, K = —1.9. The connecting dashed 
lines are for visual aid only. 



When SPF is used to control chaos, system parameter 
a is perturbed by an amount Aa proportional to the error 
Ax from the fixed value. The equations for the system's 
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FIG. 6: Data (dots) and prediction (open circles) showing the 
first difference of system values d(xi) = (xi — Xi-i) for non- 
steady sinusoidal convergence, with feedback gain K = —2.3. 
Connecting dashed lines are for visual aid only. 
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FIG. 7: Data (dots) and prediction (open circles) showing the 
first difference of system values d(xi) = (xi — Xi-i) for non- 
steady sinusoidal convergence, with feedback gain K = —3.4. 
Connecting dashed lines are for visual aid only. 



Eq. (gj except for J 2 i = Kf x . This makes det(J) = 
resulting in stability when — 1 — f x < Kf a < 1 — f x 
and solution Ax n = A(f x + Kf a ) n . Thus the behavior 
of the transient response for SPF is quite different from 
that for DFC. The SPF transient response has a single 
characteristic multiplier so it is a steady geometric con- 
vergence, either alternating sign or monotonic depending 
on whether (f x + Kf a ) is positive or negative. Immediate 
convergence occurs when K — —f x /f a since this causes 
the characteristic multilplier to be zero. DFC has either 
a more complicated geometric convergence governed by 
two characteristic multipliers, or the sinusoidal conver- 
gence that is responsible for erratic appearing non-steady 
convergence. DFC can not have immediate convergence 
because according to Eqs. (|10l) and (fT4"]) r in Eq. (fT5|) can 
never be zero. The fastest convergence occurs when K f a 
has its minimal value given by Eq. (1101) . 

For SPF applied to the Henon map with ao = 1.9 
the value of K for immediate convergence is —f x /fa — 
-1.93/0.259 = -7.45. The range of K that give stability 
is (-1 ~ fx)/ fa = -3-59 to (1 - f x )/f a = -11.3. For 
DFC, Eq. ([6]) predicts that the range of gain giving sta- 
bility is —1.8 to —3.87. This is reasonable since the error 
* ) can be expected to typically be about half of 
the first difference (x n — x n -i), so that the SPF method 
needs a gain about twice that of the DFC gain in order 
to get a similar perturbation Aa. 

We have shown that there are a variety of tran- 
sient responses when applying delayed feedback control 
to a chaotic system governed by a finite difference 1- 
dimcnsional map. In addition, the behaviors of the con- 
vergence are distinctly different from the simpler geomet- 
ric convergence obtained with simple proportional feed- 
back. Knowledge of this variety is important in situa- 
tions in which the system's response to application of 
a control algorithm is being closely monitored since er- 
ratic and non-steady convergence could be mistaken for 
a faulty control mechanism. The important advantage 
DFC has over SPF is that DFC uses an error signal that 
does not depend on the unstable fixed point. This sim- 
plifies experimental implementation and also means that 
DFC will continue to stabilize chaos even if the unsta- 
ble fixed point drifts by small amounts due to noise or 
changing system parameters. 



value and parameter are, 

x n+1 = f (x n ,a n ) (21a) 
a n+ i = a Q + K (/ (x n ,a n ) - x*) . (21b) 

The Jacobian for linear stability analysis is the same as 
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